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ABSTRACT 

Gamma-ray bursts (GRBs) are luminous enough to be detectable up to redshift z ~ 10. 
They are often proposed as complementary tools to type-la supernovae (SNe la) in 
tracing the Hubble diagram of the Universe. The distance calibrations of GRBs usually 
make use one or some of the empirical luminosity correlations, such as ri ag — L, V — L, 
E p — L, E p — Ej, 7rt — L and E p — E\ so relations. These calibrating methods are based 
on the underling assumption that the empirical luminosity correlations are universal 
over all redshift range. In this paper, we test the possible redshift dependence of 
six luminosity correlations by dividing GRBs into low-z and high-;? classes according 
to their redshift smaller or larger than 1.4. It is shown that the E p — E 7 relation 
for low-2 GRBs is consistent with that for high-z GRBs within ler uncertainty. The 
intrinsic scatter of V — L relation is too larger to make a convincing conclusion. For 
the rest four correlations, however, low-z GRBs differ from high-z GRBs at more 
than 3er confidence level. As such, we calibrate GRBs using the E p — E 1 relation in 
a model-independent way. The constraint of high-z GRBs on the ACDM model gives 
VLm = 0.302 ± 0.142(la), well consistent with the Planck 2015 results. 
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1 INTRODUCTION 


The consistent lumi nosity of type-la s upernovae (S Ne l a) ma kes them to be standard candles in probing the expansion 

1 19981 : Perlmutter et al. 19991 ). However, the redshift of SNe la is usually less than 1.4. 


history of the Universe (|Riess et al 


Gamma -ray bursts (GRBs), as the most energetic explosions in the Universe, are brigh t enough to be detectable up to redshift 
z ~ 10 ( Salvaterra et al. 20091 : Tanvir et al. 2009; C ucchiara et al. 2011 : Tanvii 20131 ). Therefore, they are often proposed as 


complementary tools to SNe la in tracing the Hubble diagram of the high-redshift Universe. Actually, GRBs have already 


(Schaefer) 2003 

Bloom, Frail fc Kulkarni 2003: Xu. Dai & Liang 200,4 Firmani et al. 2005j: Liang Rr. Zhang 2005: Firmani et al. 

200f|schaefer 

200?!: Liang et al. 20081; Liang & Zhang 200s!; Wei & Zhang 2009: Wei 201ol; Demianski & Piedinalumbol 201ll: 

Wang, Qi & Dai 2011; Capozziello et al. 

2012 

Amati & Della Valle 20l4 Wei, Wu & Melia 2013.; Velten. Montiel & Carneiro 

2013 

Cai et al. 2013 

: Breton & Montiel 

2013 

Chang et al. 2014: Cano & Jakobsson 2014: Cuzinatto. Medeiros & de Morais 

2014; 

Wang & Wang 

20141: Wang, Dai & Liang 20l4 Li, Ding & Zhu 2015). Unfortunately, since the explosion mechanism of 


Many methods have been proposed to calibrate GRBs (Dai, Liang; Sz Xu 2004i; Ghirlanc 

a et al. 2004!;; Liang; $z Zhang] 2005 

Firmani et al. 20051: Schaefer 2007: Liang et al. 20081: Liang & Zhang 20081: Wei fe Zhang 

2009; Wei 2010; Liu & Wei 2014]) 


Most calibrating methods rely on the empirical luminosity correlations found in long GRBs. At least six luminosity correlations 
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can be used in the calibration. Nom^Marani^^onnell ( 200 o|) found a correlation between spectrum lag and isotropic peak 
luminosity (ri ag — L relation). F enimore fc Ramirez-Ruiz ( 20001 ) found a correlation between time variability and isotropic 
peak luminosity (V — L relation). Amati et al. I ( 2002 ) found a tight correlation be t ween the peak energy of vF u spectrum and 
isotropic equivalent energy ( E p — Ei so relation). Ghirlanda. Chisel ini fc Lazzati (12004 1 found a similar correlation between 
peak energy and collimation-corrected energy (E p — E-, relation), ^netoku et all (|2004 1 found a correlation between peak 
energy and isotropic peak luminosity ( E p — L relation). Schaefer 1 2007{ l found a correlation between minimum rise time of 
light curve and isotropic peak luminosity (trt — L relation). 

All of the calibrating methods based on the empirical luminosity correlations have an underlying assumption, that is, 
the luminosity correlations do not evolve with redshift. If the luminosity correlations is not universal over the whole redshift 
range, these calibrating met hods w ill fail. In fact, the possible redshift dependence of luminosity correlations has already been 
tested by some authors. Basilakos fc Perivolaropoulosl ( 2008 1 investigated the above six empirical luminosity correlations in 
four redshift bins, and showed that the slopes of all six correlations differs between redshift bins, although the intercepts do 
not vary significantly. Since the CRB sample is not large enough in each bin, the statistical uncertainty is large. Therefore, 
they concluded that no statis ticall y signi f icant evidence for the redshift evolution of the luminosity correlations was found. 
With the updated data, Wang, Qi fe Dai (2011) got a similar conclusion. However. [Li (2007) investigated the Amati relation 
in four redshift bins and showed that the slope and intercept varies with redshift systematically and significantly. Recently, 
Lin et alJ (201a) divided GRBs into two redshift bins, and found that the Amati relation (espec ially the slope parameter) 
of low -0 GRBs differs from that of high-z GRBs at more than 3cr confidence level. Dainotti et al. (2013;) investigated the 
slope evolution of GRB correlations and showed that correlation slope that differs from the intrinsic one may overestimate or 
underestimate the cosmological parameters. 

In this paper, we recheck the possible redshift dependence of six luminosity correlations. We divide GRBs into low-z and 
high-z classes according to their redshift smaller or larger than 1. 4, and test the lum inos ity correlations for low-z and high-z 
GRBs, respectively. The main difference between our work and Wang. Qi fc 53 (l 201 lh ’s is that we just divide GRBs into 
two redshift bins, so that the number of GRBs in each bin is large enough to do statistical analysis. We choose z = 1.4 as 
the threshold because the redshift of SNe la is usually smaller than 1.4, and the Universe below this redshift has already 
been tightly constrained. We find that, among the six luminosity correlations, only the E p — E 7 relation is consistent between 
low-z and high-z GRBs within l a unce rtainty. As such, we can calibrate GRBs through the E p — E 1 relation using the Pade 
approximation proposed by LhjL& Wei ( 2014T L and the Hubble diagram of GRBs can be constructed. 

The rest of the paper is arranged as follows: In section [2] we test the redshift dependence of six luminosity correlations. 
In section [3] we calibrate the distance of high-z GRBs using the E p — E 1 relation, and then use them to constrain the ACDM 
model. Finally, a short summary is given in section [4] 


2 TESTING THE REDSHIFT DEPENDENCE OF LUMINOSITY CORRELATIONS 


All the six luminosity correlations mentioned above have the exponential form R = AQ b , which can be linearized by 
taking the logarithm, i.e., 

y = a + bx (y=logR, a; = log Q, a = log A), (1) 

where “fog” represents the logarithm of base 10. For the sake of clarity, we write the six luminosity correlations explicitly 
here: 
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where quantities with a subscript “i 

\-i 


represent the quantities in the comoving frame, which can be transformed to the 
observer frame by Ti agii = ri ag (l + z) -1 , r R T,i = tr T (1 + z) -1 , V) = V(1 + z), and E Pti = E p ( 1 + z). 


The isotropic peak luminosity L can be calculated from the bolometric peak flux Pboio as Schaefer! 20071 ) 
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where di is the luminosity distance. The bolometric peak flux Pbolo is c alculated from th e observed peak photon flux in the 
rest frame 1 — 10, 000 keV energy band by assuming the Band spectrum ( Band et al. 19931 ). The luminosity distance depends 
on a specific cosmological model. In the concordance ACDM model, it is given as 

dz 


d L (z) = (l+z) 


Ho 


F 

JO 


(9) 


\ZOm(1 + z)^ + (1 — Dm) 

where Dm is the mater density, Ho is the Hubble con stant, a nd c is the light speed. Here we take Dm = 0.280 and Ho = 
70.0 km s -1 Mpc -1 from fitting to the Union2.1 dataset ( Lin et al .203). The uncertainty of L propagates from the uncertainty 
of Pbolo, while that from db is a bsorbe d i nto t he intrinsic scatter. The isotropic equivalent energy Ei so can be calculated from 
the bolometric fluence Sholo as ([S chaefer 1200711 




= 47 rd| 5 boio(l + z) \ 


( 10 ) 


Similar to the bolometric peak flux, the bolometric fluence Sbolo also corresponds to the rest frame 1 — 10, 000 keV energy 
band. For Ei so , we also only consider the error propagation from Sbolo. The collirnation-corrected energy, P 7 , is the isotropic 
equivalent energy multiplied by a beaming factor Pbeam = 1 — cos 9jet, where 9j e t is the jet opening angle, i.e, 

E~f = PisoPbeam = 47rd|SboloPbeam (1 + z) _1 - (11) 


The uncertainty of P 7 propagates from the uncertainties of both Sbolo and Pbeam- The error propagation from Q to logQ is 
given as 


1 OQ 

aiosQ ~ MIOQ - ’ 


( 12 ) 


where “In” represents the natural logarithm. If Q has nonsymmetric error, we symmetrize it by taking the average, i.e., 
a Q = (°q + *q)/2. _ 


To test the possible redshift dependence of luminosity correlations, we analyze the GRB sample taken from Wang . Qi fc Dai 


(2011J). This sample consists of 116 long GRBs in the redshift range 2 £ [0.17,8.2]. This dataset is a collection of GRBs with 
well-measured spectra properties from various instruments, such as BATSE, Konus, Swift, etc.. We divide GRBs into two 
subsamples according to their redshift smaller or larger than 1.4, and call them low -2 and high -2 subsamples, respectively. We 
choose 2 = 1.4 as the threshold because th e redshift of SNe la is usually smaller than 1.4. The Universe below this redshift 
has already been well studied using SNe la ( Amanullah et al. 20101 : Suzuki et al. 2012 : Betoule et al. 2014 ). Low -2 and high -2 
subsamples consist of 50 and 66 GRBs, respectively. We fit each luminosity correlation to the two subsamples separately. 
Since the plot of each correlation in the xy plane show large error bars in both the horizontal and vertical axes, and intrinsic 
scatter domi nates over t h e me asurement error, the ordinary least-X ' 2 method does not work well. We apply the fitting method 
presented in D’Agostini ([200 51). The best-fit parameters (a,b,a i n t) can be derived by maximizing the D’Agostini’s likelihood, 


Pu(o-int,a, b) oc 


V° 


+ b 2 c 


x exp 


(yi — a — bit ) 2 


20 


°Tnt + 


+ b 2 a%.) 


(13) 


where the intrinsic scatter dint represents any other unknown errors except for the measurement error. Equivalently, we can 
minimizing the x 2 , 


Xd(^ int, a,b) = ^2 ln(<r 2 nt + a 2 , + & 2 a 2 .) + ^ 


(yi - a - bxj) 2 
ofnt + o/L + b 2 a 2 


(14) 


We use the publicly available Matlab package FMINUIT0 to derive the best-fit parameters and their uncertainties. The 
results are listed in the fourth to sixth columns in Table [T] This table gives the mean values of the best-fit parameters and 
their la uncertainties. Note that not all GRBs are available in the analysis of each luminosity correlation. For example, GRBs 
without measurement of the jet opening angle is unavailable in the E p — P 7 analysis, while GRBs having no spectrum lag 
measurement are invalid in the ri ag — L analysis. For this reason, we also list the number of available GRBs in each fitting in 
the third column of Table [T] All the six luminosity correlations are plotted in Figure [T| in logarithmic coordinates. Low -2 and 
high -2 GRBs are denoted by black and red dots, respectively. The error bars represent lcr uncertainties. Since the Swift/BAT 
instrument is only sensitive in a narrow energy band (~ 15 — 150 keV), the uncertainties of peak energy of some GRBs are 
extremely large. The lines stand for the best-fit results (black line for low -2 GRBs and red line for high -2 GRBs). Besides, we 
also plot the lcr, 2cr and 3 a contours in the (a, b) plane for low -2 (black curves) and high -2 (red curves) GRBs in Figure [2] 
The best-fit central values are denoted by dots. 

From Tableland Figure[T] we can see that among six luminosity correlations, the V — L relation has the largest intrinsic 
scatter, while the E p — P 7 relation has the smallest intrinsic scatter. The intrinsic scatter of the V — L relation is so large that 
it is unreasonable to fit it with a line. For all the six luminosity correlations, high -2 GRBs have larger intercept, but smaller 
absolute slope than low -2 GRBs, although the difference of intercepts between low -2 and high -2 GRBs is not as significant 


1 http: //www. fis.unipr.it/~giuseppe.allodi/Fminuit/Fminuit_intro.html 
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Table 1. The intrinsic scatters (ai n t), intercepts (a) and slopes (6) of six luminosity correlations for low -z and high -2 GRBs, derived 
from maximizing the D’Agostini’s likelihood. The quoted errors are of I a. N is the number of GRBs available in the fitting. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

correlation 

subsample 

N 

Tint 

a 

b 

"Hag 

low-z 

27 

0.475 ± 0.066 

52.103 ± 0.093 

-0.783 ±0.141 

high-z 

32 

0.321 ±0.049 

52.486 ± 0.064 

-0.633 ± 0.111 

V-L: 

low-z 

47 

0.875 ±0.101 

51.550 ± 0.179 

0.441 ± 0.266 


high-z 

57 

0.546 ± 0.057 

52.262 ± 0.121 

0.221 ±0.117 

Ep — L\ 

low-z 

50 

0.577 ±0.063 

51.880 ± 0.087 

1.461 ±0.174 


high-z 

66 

0.386 ± 0.039 

52.363 ± 0.054 

1.102 ±0.143 

Ep — E 1 : 

low-z 

12 

0.159 ±0.059 

50.637 ±0.061 

1.539 ±0.144 


high-z 

12 

0.261 ±0.105 

50.649 ± 0.096 

1.354 ±0.275 

trt — L: 

low-z 

39 

0.470 ± 0.058 

52.685 ± 0.122 

-1.318 ± 0.181 


high-z 

40 

0.395 ± 0.052 

52.747 ± 0.078 

-0.784 ±0.158 

Ep E[ so : 

low-z 

40 

0.561 ±0.069 

52.561 ± 0.095 

1.586 ±0.190 


high-z 

61 

0.365 ± 0.040 

52.874 ± 0.053 

1.243 ±0.134 


as that of slopes. The slope difference of the trt — L relation is especially evident. This can be seen more clearly from the 
contour plots in the (a, b ) plane in Figure [2] The E p — E 7 relation of low-z GRBs is consistent with that of high-z GRBs 
within la uncertainty. However, for the rest five luminosity correlations, low-z GRBs differ from high-z GRBs at more than 
3cr confidence level. Especially, there is no overla p betw een the 3a contours of two subsamples for the E p — L relation. As for 


the Amati relation, we recover the results of Lin et al. 1 2015h . 


The results above are derived using D’Agostini’s likelihood. Since the observed data points. ha ve sig nificant errors on 
both the x-axis and y-axis, there is no unique method to determine the best-fit parameters. Reichartl ( 211(11 ) has constructed 
a likelihood which is slightly different from D’Agostini’s one. To test whether the above results depend on the choice of a 
specific best-fit method, we also do a similar calculation using Reichart’s likelihood. The Reichart’s likelihood is written as 
( Reichart 20011 


£r(g 


±* n 


CTTTZ 


Ctf + + b 2 ( CT I + a Zi ) 


x exp 


{yi-a - bxi) 2 


2[a2 + a 2 


+ b 2 {a 2 ± erg.)] 


where a x and a y are the intrinsic scatters along the x-axis and y-axis, respectively. The corresponding \ 2 is given as 


Xr(<7x 


,a,b) = 5^1n[o 


+ b 2 (al ± al.)] + ^2 


(yi-a - bxi) 2 


a y + a m + b 2 {<r% + cr|.) 


- Nln{l + b 2 ), 


(15) 


(16) 


where N is the number of data points. The best-fit parameters are the one which can minimize the right-hand-side of Eq. m- 

The best-fit parameters and their lo- uncertainties are listed in Table [2] The last column gives the “equivalent” total 
intrinsic scatter, which is calculated from o-i n t = (a^ + Cal) 1 ^ 2 . Comparing to Table]!] we can see from Table[2]that the E p — E 1 
relation has the smallest (while the V — L relation has the largest) intrinsic scatters, although the uncertainties of intrinsic 
scatters in Table [2] are much larger. Using Reichart’s likelihood, the parameters (especially the intrinsic scatter) cannot be 
well constrained. Reichart’s likelihood leads to larger absolute slope parameters compared to D’Agostini’s likelihood. Figured] 
is the contour plot in the (a, b) plane. We can see an important common feature between the results derived from two different 
likelihoods: only for the E p — E 1 relation, low-z subsample is consistent with high-z subsample within lcr uncertainty. For 
the Ti ag — L and E p — L relations, low-z subsample still differs from high-z subsample at more than 3cr confidence level. As 
for the trt — L and E p — Ei SO relations, low-z subsample differs from high-z subsample at more than 2a confidence level. The 
relatively lower significance is due to the larger uncertainties of the best-fit parameters. The uncertainties of slope parameters 
of V — L relation derived from Reichart’s likelihood are extremely large. In a word, only E p — E~, relation shows no significant 
evidence for the redshift evolution. This conclusion does not depend on the choice of the best-fit methods. 

Reichart’s likelihood differs from D’Agostini ’ s one by an extra factor (1 + b 2 ) 1 ^ 2 . Otherwise, these two likelihoods are 
identical (if we set a ; 2 nt = a 2 + b 2 al). D’Agostinfl (2005) pointed out that Reichart’s likelihood has a problem: m 2 cannot be 
added tout court to 1, since m 2 is in general dimensional (although in our case it is dimensionless). The factor (1 + b - 2 ) 1 / 2 has 
the net effect of overestimating m. This is one reason why Reichart’s likelihood leads to a larger slope parameters relative to 
D’Agostini’s likelihood. Therefore, we use the results of D’Agostini’s likelihood when calibrating the distance of GRBs in the 
next section. 
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log(-B p ,i/300 keV) 





log(-Ej,,i/300 keV) 


Figure 1. The luminosity correlations for low -2 (black) and high -2 (red) GRBs. Error bars represent the la uncertainties. The lines are 
the best-fit results, which are derived from maximizing the D’Agostini’s likelihood. 


3 DISTANCE CALIBRATION AND COSMOLOGICAL IMPLICATIONS 


As we have shown that the E p — U 7 relation does not signi ficantly evo l ve wi th redshift, we can use it to calibrate GRBs. 
To avoid the circularity problem, the Pade method proposed by Liu fe Weil ( 2014 ) is applied. The main cali bra ting p rocedures 
are as follows: Firstly, derive the distance-redshift relation of SNe la (here we use the Union2.1 ( Suzuki et al. 2012 ) dataset) 
using the Pade approximation of order (3,2), i.e., 


a{z) = 


ao + OL\Z + Ct2Z 2 + «32 3 

l + /3i z + p 2 z 2 


(17) 

















































































6 H.-N. Lin, X. Li and Z. Chang 



Rag L 



E p — L 



Trt — L 



V-L 

2.5 
2 

1.5 
1 

0.5 
0 

50.2 50.4 50.6 50.8 51 



E p — Ej, 



Figure 2. The la, 2<r and 3<x contours in the (a, b ) plane for low -2 (black curves) and high-z (red curves) GRBs derived from the 
D’Agostini’s likelihood. The central values are denoted by dots. 


where the coefficients (ap, ct i, Q2 , <33 , /3i , ft) and the corresponding covariance matrix are derived by fitting Eg. (11711 to the 
Union2.1 dataset (see Liu fe We] ( 2014 ) for details). Assuming that the low-z GRBs trace the same Hubble diagram to SNe 
la, we can calculate the distance moduli of low-z GRBs directly from Eq. EJ- The uncertainty of /r propagates from the 
uncertainties of the coefficients (a;, ft). Then the luminosity distance of low-z GRBs can be obtained using the relation 


H(z) = 5 log + 25. (18) 

As d,L is known, the collimation-corrected energy can be further calculated from Eq. CEB- Note that there are only 12 low-z 
GRBs and 12 high-z GRBs available since the others have no measurement of jet opening angle. Then we fit the E p — E-y 



















Distance calibration of GRBs and constrain on the ACDM 7 


Table 2. The intrinsic scatters along the rr-axis (<r x ) and y-axis (<Jy), intercepts (a) and slopes (b) of six luminosity correlations for low- 2 : 
and high- 2 : GRBs, derived from maximizing the Reichart’s likelihood. The quoted errors are of lcr. N is the number of GRBs available 
in the fitting. The last column gives the total intrinsic scatter cri n t = (a 2 + b 2 ^ 2 ) 1 ' 2 . 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

correlation 

subsample 

N 

&X 

Oy 

a 

b 

°mt 

Rag L: 

low -2 

27 

0.271 ± 1.757 

0.417 ± 1.385 

52.136 ±0.102 

-1.101 ± 0.183 

0.513 ± 1.593 


high- 2 

32 

0.236 ± 1.390 

0.255 ±0.855 

52.471 ±0.065 

-0.815 ± 0.115 

0.319 ±0.966 

V — L: 

low -2 

47 

0.395 ± 0.055 

0.010 ± 29.430 

47.945 ± 2.413 

6.135 ±3.571 

2.425 ± 1.457 


high-2 

57 

0.390 ± 0.754 

0.236 ±8.803 

49.730 ± 1.324 

2.661 ± 1.206 

1.064 ±2.804 

E p — L: 

low -2 

50 

0.042 ± 3.396 

0.666 ± 1.053 

51.950 ±0.102 

2.235 ±0.252 

0.672 ± 1.479 


high-2 

66 

0.091 ± 1.036 

0.396 ± 0.829 

52.309 ± 0.065 

1.871 ±0.208 

0.431 ± 1.078 

Ep — E'yl 

low -2 

12 

0.047 ± 1.121 

0.152 ±0.951 

50.634 ± 0.064 

1.663 ±0.163 

0.171 ± 1.198 


high-2 

12 

0.159 ± 0.931 

0.073 ± 6.867 

50.598 ±0.119 

1.836 ±0.399 

0.300 ± 2.347 

Trt — L: 

low -2 

39 

0.073 ± 1.438 

0.551 ±0.880 

53.088 ±0.182 

-2.144 ±0.317 

0.573 ± 1.196 


high-2 

40 

0.222 ± 0.778 

0.293 ± 1.291 

52.873 ± 0.094 

-1.481 ± 0.270 

0.440 ± 1.217 

E p Ei so : 

low -2 

40 

0.225 ± 0.658 

0.373 ± 2.042 

52.631 ±0.108 

2.270 ±0.255 

0.632 ± 1.705 


high-2 

61 

0.160 ±0.474 

0.275 ± 0.955 

52.838 ± 0.063 

1.863 ±0.193 

0.405 ±0.917 


Table 3. The distance moduli of 12 high -2 GRBs calibrated through the E p — E 7 relation. 


GRB 2 


010222 

1.48 

030328 

1.52 

990123 

1.61 

990510 

1.62 

030226 

1.98 

021004 

2.32 

050820A 

2.61 

030429 

2.66 

050401 

2.9 

020124 

3.2 

060526 

3.21 

060605 

3.8 


A 1 


45.0990 

0.4684 

44.8357 

0.5404 

45.5094 

0.6244 

45.6337 

0.4921 

46.6219 

0.6476 

46.0496 

0.9760 

47.4273 

0.7877 

46.8414 

0.8608 

47.1522 

0.6272 

46.5993 

0.6032 

45.4623 

0.6342 

50.9127 

1.0548 


relation (i.e., Eq.©) to the 12 low -2 GRBs, which gives the best-fit parameters 

a int = 0.161 ±0.059, a = 50.632 ± 0.062, b = 1.537 ± 0.145. 


(19) 


By directly extrapolating the E p — E 7 relation to high -2 GRBs, we can inversely obtain the collimation-corrected energy for 
12 high -2 GRBs from Eq. J5]). Finally, calculate the luminosity distance of high -2 GRBs from Eq. CD, and then the distance 
moduli fr om Eq. juj). The uncertainty of distance moduli propagates from the uncertainties of E 7 , Sboio and fb ea m, i.e. 
( Schaefer 2007 1. 


2 In 10 


(In 10)"cri 2 O gB 7 + 


'-’bolo 

02 

°bolo 


+ 


^^beam 


F 2 

beam 


where 


2 

^log E-J 


= CT a + l a b log 


E v 


300 keV 




In 10 E v 


I 2 

± (Rnt • 


( 20 ) 


( 21 ) 


The distance moduli of 12 high -2 GRBs and their lcr uncertainties calibrated through the E p — E 1 relation are listed in 
Table [3] We also plot the 12 high -2 GRBs in z — fi plane in Figure [4] where the black curve is the best-fit result to ACDM 
model. The fit of 12 high -2 GRBs to the ACDM model gives !1m = 0.302 ±0.142, well consistent with the Planck 2015 results 
( Ade et al. 2015l l . From Figure 0] we can see that the distance of GRB 060605 is much overestimated. This is because the 
E p — E 1 relation overestimates the energy of GRB 060605 (see also the E v — E-y plot in Figure|T| where the red star represents 
this burst). On the contrary, the distance of 060526 is underestimated because the E p — E~, relation underestimates its energy. 
The rest 10 GRBs are consistent with the ACDM model within lcr uncertainties. 

For comparison, we also calibrate GRBs through the E p — Ei so relation (the so called Amati relation). In this case, 40 low -2 
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Figure 3. The la, 2a and 3c contours in the (a, b) plane for low -2 (black curves) and high-z (red curves) GRBs derived from the 
Reichart’s likelihood. The central values are denoted by dots. 


GRBs and 61 high-z GRBs are available. The constraint of 61 high-z GRBs on the ACDM model gives 12 m = 0.805 ± 0.144, 
which is much larger than the Planck 2015 results. The reason for this can be easily understood. From the E p — Ei so plot in 
Figure [H we can see that high-z GRBs have in average larger isotropic equivalent energy than low- 2 : GRBs at the same E p 
value. Therefore, when extrapolating the Ainati relation from low -2 GRBs to high -2 GRBs, the energy (so the distance) of 
most high -2 GRBs is underestimated. The underestimation of distance further leads to the overestimation of 12 m. For this 
reason, we can predict that GRBs calibrated through the rest four luminosity correlations (n ag — L, V — L, E p — L and 
trt — L) may also overestimate the value of 12 m • 
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Figure 4. The Hubble diagram of 12 long GRBs calibrated using the E p — E-y relation. The black curve is the best-fit result to ACDM 
model. The best-fit parameter is Qm = 0.302 ± 0.142. 


4 SUMMARY 


In this paper, we checked the possible redshift dependence of six luminosity correlation in long GRBs. We divided GRBs 
into low -2 and high -2 subsamples according to their redshift smaller or larger than 1.4. The slope and intercept parameters 
of six luminosity correlations are derived by maximizing the D’Agostini’s likelihood. For all the six luminosity correlations, 
high -2 GRBs seem to have larger intercept, but smaller absolute slope than low -2 GRBs. It was shown that the intrinsic 
scatter of V — L relation is to large to make a convincing conclusion. The E p — E~, relation has the smallest intrinsic scatter 
among the six, although the number of available GRBs is small. Most importantly, the E p — E~, relation shows weak redshift 
dependence. Strong evidence (> 3 <j) for the redshift evolution was found in the rest four correlations. Similar features can 
be seen when we use Reichart’s likelihood instead of D’Agostini’s, although the statistical significance is lower. We calibrated 
high -2 GRBs using the E p — E 1 relation in a model independent way and reconstruct the Hubble diagram. The constraint of 
high -2 GRBs on the ACDM model gives matter density SIm = 0.302 ± 0.142, which is well consi stent with the Planck 2015 
results, although the error bar is large. Calibrating GRBs using the Amati relation, as was done by Liu & Wei (|2014l ). in some 
cases may overestimate Dm- One of the disadvantage in using the E p — E 1 relation, of course, is that only a small number 
of GRBs are available since most GRBs have no measurement of jet opening angle. We hope that the future observation will 
enlarge the GRB sample so as to improve the statistical significance. 
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